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ABSTRACT 



Aims. We propose a method for measuring the projected rotational velocity vsini with high precision even in spectra with blended 
lines. Though not automatic, our method is designed to be applied systematically to large numbers of objects without excessive com- 
putational requirement. 

Methods. We calculated the cross correlation function (CCF) of the object spectrum against a zero-rotation template and used the 
Fourier transform (FT) of the CCF central maximum to measure the parameter v sin i taking the limb darkening effect and its wave- 
length dependence into account . The procedure also improves the definition of the CCF base line, resulting in errors related to the 
continuum position under 1 % even for v sin i = 280 km s~ 1 . Tests with high-resolution spectra of F-type stars indicate that an accuracy 
well below 1% can be attained even for spectra where most lines are blended. 

Results. We have applied the method to measuring vsiin in 251 A-type stars. For stars with vsin; over 30 kms~' (2-3 times our 
spectra resolution), our measurement errors are below 2.5% with a typical value of 1%. We compare our results with Royer et al. 
(2002a) using 155 stars in common, finding systematic differences of about 5% for rapidly rotating stars. 

Key words, stellar rotation; astronomical techniques; A-type stars 



1. Introduction 

The projected axial rotational velocity of single stars can be mea- 
sured directly from the broadening of their spectral lines. That is 
why the rotational velocity is an important observable in statis- 
tical studies of stellar astrophysics. Several methods have been 
developed to measure v sin z, but the problem of systematic dif- 
ferences between different authors and methods has always been 
present. Th e first workable mode l for the stellar rotation was es - 
tablis hed bvlCarroll I dl928l[T933b and lCarroll & Ingram I (fl933l) . 
but Shai n & Struve I (Il929h presented a very simple graphical 
model called the classical model of a rotating star (CMRS) by 
y— ( , ICollins & Truax I |l995), which was used as a standard for most 
I* ■ of the work done in the 2 0th century (for more details see the ex- 
. — cellent review paper by ICollinsl |2004|) . IShain & Struve l (fl929T) 
. did not co nsider limb dar kening, which was introduced in the 
r Cl 1 CMRS bv lCarrolll(ll933h . In a series of papers Sletteba k mea- 
?3 \ sured v sin i for stars in the spectral rang e O-G (ISlettebaklll95i 
ll955Ul956tlSlettebak&Howardlll955l) by considering a linear 
limb darkening law. These stars were used in many subsequent 
papers by different authors to calibrate the full width at half 
depth of the line used to measure the v sin i parameter. Therefore, 
determination of v sin i was conditioned by the calibration made 
by the author. 

In the last quarter of the 20th centurv lsiettebak et alTI dl975l) 
established a new standard system of stars distributed over both 
hemispheres and covering a range in spectral type from 09 to F9. 
The measure of v sin i was again based on a calibration of v sin i 
with the full wi dth at half depth of st ellar absorption lines in the 
star spectrum. Slettebak et al. I d 1975b established this system by 
comparing their digital data with numerical models constructed 



by adopting Roche geometty, uniform angular velocity at the 
surface, von Zeipel gravity darkening, and numerical integration 
of angle-dependent model atmosphere intensities. The standards 
for th e high-velocity specimens were established by Slettebak 
(1982) by visual comparison of line widths on the spectrograms 
that were recorded on photographic plates. Systematic effects 
between the old and ne w Slettebak systems were studied by 
iGarcfa & Levato I ((1984). The calibrations of the systems de- 
scribed were made for only three lines, He 1 4471 A, Mgn 4481 
A, and Fe i 4476 A, and the theoretical profiles were calculated 
with main sequence models. Besides this, determination of the 
full width at half depth has shown to be very sensitive to the 
continuum position, which represents an important error source. 

A significant improvement was possible thanks to methods 
based on Fourier transform (FT) of line profiles providing a 
deeper analysis. In the first place, this tool is used to avoid an 
external calibration by suppressing the error related to this stage 
of the process. In second place, with high signal-to-noise ra- 
tio (S/N) it is possible to identify another broadening agents 
present in the line profile, at least in a qualitative way. Finally, 
it is possi ble to ana l yze second-order effects like differential 
rotation (lGravlll977lll982l:lBmning1ll98ll:lReiners & Schirnttl 
2002 ) and to put a Umit on the inclination of the rotational axes 
dReiners & Rover1 l2004). Regardless of the method adopted to 
determine v sin i, it is not an easy task to find the spectral lines 
that have the minimum conditions required to be used in the 
measurement of rotational velocity, namely: 1) lines with no 
blends, 2) lines intense enough to be identified in rapid rota- 
tors, and 3) lines mainly broadened by rotation. For instance, in 
O-type and B-type stars, only the Balmer lines and some helium 
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lines are intense enough, but they all show a significant Stark 
effect. In stars of spectral type later than A4, almost all lines 
are blended if vsin/ is grater than 100 kms" 1 , making it im- 
possible to_jmd_anj£olatejHine. This problem is evident in the 
work of iRover et al. I d200 2a) where fewer than three lines were 
measured in A-type stars with vsin / > 60 kms -1 . Because this 
problem increases with the spectral type, other methods have 
been developed. 



One solution to the blending limitation is the use of a least 
square deconvolution procedure to derive the broadening func- 
tion in a selected wavelength region instead of a single line pro- 
file. This methodology implies the application of an iterative 
process to fit the equivalent width of the template's spectral lines, 
the broadening function, and the continuum position (for details 
see lReiners & Schmitt 1 12003). This is a very powerful technique 
for a detailed study of the rotational profile. However, for exten- 
sive applications like the development of a catalog of rotational 
velocities, a more direct method that does not involve fitting the 
intrinsic spectrum or including any atmospheric parameter other 
than limb darkening, might be more suitable. 

Even though CCF had been originally proposed for deter- 
mining radial velocities, the projected axial rotational velocity 
can be inferred from the width of the CCF maximum. This tool 
has often been applied to determine v sin i in cool stars (M-L 
spectral type). The standard procedure is to calculate the CCF 
between an observed spectrum and a template spectrum of the 
same temperature with vsin / = kms 1 . Then, a fitting profile 
for the maximum of this function is calculated and the width of 
the fit can be used to measure v sin / through an empirical calibra- 
tion v sin /-width. The absence of single lines is somehow solved 
by means of the CCF . Nevertheless, as in any other method that 
depends on empirical calibrations, various strategies, not always 
equivalent, have been adopted in the literature. I n some works 
the c entral maximum is fitted with a Gaussian dBailer- Jones I 
|2004|) . while in others a parabola dTinne v & Reid Il99q). or even 
a Gau ssian plus a quadratic function are used dWhite & Basri 
120031) . The adoption of different functions to evaluate the rota- 
tional broadening could lead to systematic differences in results 
from different authors. 

Templates selection is also very heterogeneous. Even though 
some authors calculate synthetic template spectra, real stellar 
spectra have been used in most works based on the CCF (see 
Whit e & Basril 120031; iMohantv & Basri I 120031: iBailer-JonesI 
2004 . Since the CCF contains information from the template 
and the object spectrum, using an observed template could in- 
troduce an external error source from the unknown broadening 
factors present in the template spectral lines. Nowadays, the best 
methods dealing with stars that show intense line blending in 
their spectra require significant computational resources. 

Motivated by the need for a precise and expeditious tech- 
nique to be applied extensively for the construction of a catalog 
of rotational velocities of bright A-type stars (Levato et al. in 
preparation), we develop here an alternative method based on 
the CCF and, at the same time, independent of any external cal- 
ibration. Our procedure uses the FT to measure the parameter 
v sin /, taking the dependence of the transform with limb darken- 
ing into account. In Sect. 2 a full description of the methodology 
is presented. Limb darkening consideration and other practical 
details of the procedure are explained in Sect. 3. Sect. 4 describes 
specific application to A and later-type stars. The precision of 
the obtained rotational velocities is discussed in Sect. 5, and our 
main conclusions are summarized in Sect. 6. 



2. Method 

As a first approximation, according to the c l assica l model of 
a rotating star revised by [Collins & Truax I d 19951) (see also 
iBrown & Verschueren II 19971) . we can consider an observed stel- 
lar line profile as the convolution of an intrinsic line profile and 
a rotational broadening function. Then, the observed spectrum 
D(A) can be approximated by a convolution between a template 
spectrum T(A) and the rotational broadening function G{A), 

D(A) = T(A) * G(A), 

where T(A) includes any other broadening effect different from 
rotation: natural line broadening, thermal broadening, microtur- 
bulence, Stark effect, etc. 

Under the previous assumption the CCF between the object 
spectrum and a template of the same spectral type, but without 
rotation, results in 



CCFnr - 



D*T 



(T*G)*T 



D\ 



- CCF 'tt * G, 



(1) 



where || T || and || D || are the norm of each spectrum. Since 
the rotational velocity is the main broadening agent of metallic 
lines in stars wit h vsin / > 10km s -1 , i. e . in most normal main- 
sequence stars (lAbt&Morrem Il995t IXbT Levato & Grosso 
120021) . the intrinsic width of a line profile <x is usually <x <k 
v sin i. 

The autocorrelation function CCF Yr presents a narrow peak 
centered at zero, which for the current application can be con- 
sidered as a Gaussian whose width is V2 times the width of the 
template lines. Side lobes are present in this function, but they 
are usually very small in comparison with the central peak since 
their intensity is inversely proportional to the number of spectral 
lines involved in the cross-correlation. Therefore, regardless the 
shape of the line profile of the template, the function CCFdt is 
very similar to the rotation function G and can be used to derive 
the parameter v sin i. 

The idea of the proposed method is to calculate the FT of the 
central maximum of the CCFot to derive v sin / from the posi- 
tion of its first zero. In fact, CCFdt is essentially the function G 
convolved with a narrow profile that has no significant impact on 
the position of the first zero of the FT. The central maximum of 
the CCF has the same information about the rotational velocity 
as a single line profile. The main benefit of utilizing the CCF is 
that having to select lines without blending, which represents a 
problem in stars later than late-A early-F, is avoided. 

The CCF peak, however, might be blended with side lobes, 
but the more lines used to calculate the CCF, the bigger the inten- 
sity difference between the central maximum and the secondary 
maxima. In the proposed method, we include a procedure to re- 
move the side lobes contribution as described in Sect. 12.11 

Finally, the S/N of the CCF peak is much higher than that of 
a single line. Therefore, using the CCF calculated from a large 
spectral region instead of individual line profiles presents ad- 
vantages over both line blending and S/N. Nevertheless, there 
are some restrictions on the size of the spectral region to be cor- 
related owing to the variation in the limb darkening coefficient 
with wavelength, which will be discussed in the next section. 

2.1. Extraction of the rotational profile 

The CCFdt would be equal to the rotational profile G only in 
the ideal case in which CCFtt has a single peak of negligible 
width (a Dirac 6 function). In practice, however, besides the fi- 
nite width of the central peak, the CCFdt contains several small 
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subsidiary peaks coming from the fortuitous coincidence of each 
spectral line of the object spectrum with different spectral lines 
in the template as the former is shifted with respect to the latter. 
These small peaks can be blended with the central peak due to 
the rotational broadening of the object, affecting the determina- 
tion of the base line of the CCF central maximum and hence the 
measurement of the rotational velocity. 

To overcome this problem, a key stage in the process is the 
remotion of secondary maxima. To this aim we consider the 
function CCF Vr as the sum of two functions: CCFItt, repre- 
senting only the central peak, and CCF2jt, which includes all 
the subsidiary peaks. Then, from eq.Q] 

CCF DT = CCF TT * G = CCFl TT * G + CCF2 TT * G. 

To calculate the contribution of the subsidiary peaks, we remove 
the central maximum from the CCF Vr to obtain CCFIjt, which 
is convolved with a first approximation of the rotational broaden- 
ing profile G\ extracted from the central peak of the CCFdt- The 
result is then subtracted from the CCFdt, thereby improving the 
determination of the base of the rotational profile. Considering 
the central peak of the CCF 1 tt to be much narrower than G, we 
find from the previous equation 



G ~ CCF DT - CCF2 TT * d- 



(2) 



Before applying the subtraction, two corrections are required: 
a radial velocity correction for both CCFs to be centered, and 
a scale correction to account for an eventual global line inten- 
sity difference between the object and the template spectra. As 
illustration, FigureQ] shows the application of this process to a 
spectrum of the A2-type star HR 892 with a synthetic template 
for T = 9000 K. 

Finally, once the base line of the central maximum has been 
improved, we fit the background in the surrounding region of 
this peak and subtract the fit to have the base of the rotational 
profile at zero. Through this process the rotational broadening 
profile G is obtained. 

In principle, reliable and precise rotational profile can be re- 
constructed from the CCF, as long as the rotational broadening is 
not comparable to other photospheric or instrumental broaden- 
ing effects. In Fig. [2] we compare the theoretical rotational pro- 
file of eq.[3]with the profile recovered from the CCF of a 100 A 
long spectrum (5350-5450 A). As shown in the lower panel, 
the only noticeable difference is the smoothing of the sharp cut 
at the edges, owing to the intrinsic spectral line width (0.08- 
0. 16 A). The example corresponds to an atmosphere model with 
T e ff = 8000 K and log g=4.0, convolved with a rotational profile 
of vsin/ = 60 km s 1 and s = 0.6. 

2.2. Calculation of v sin i 

The rotational profile of a spectral line centered on a wavelength 
A Q for a spherical star rotating as a rigid body and whose limb 
darkening law is linear with a coefficient s, is 



G(x) 



2(l-e)(l-;c 2 ) I/2 + f(l-x 2 ) 
n(l - f) 



(3) 



forx < 1, andG(x) = Oforx > 1, where jc = ln(/l//l o )-c/(vsin0- 
Expressing the rotational profile in terms of A In A instead of AA 
has the advantage of being indepe ndent of the line wavele ngth 
(cf. w ith eq. 17.12 in iGrav I [19921: eq.4 in iReiners & Schmittl 
120021) . In addition, the Doppler formula in logarithmic scale, 



4600 
A (A) 



4620 4640 




-1000 1000 2000 

RV (km/s) 

Fig. 1. Calculation of the rotational profile using cross- 
correlations. Upper panel: Object spectrum and template spec- 
trum. Lower panel, from top to bottom: a) Cross-correlation 
function template-template CCFtt (thin line) and the same 
function after removing the main peak and convolving with the 
provisional rotational profile, i.e. CCFItt * G (thicker line), 
b) The same function CCFItt * G scaled to the object inten- 
sity, c) Cross-correlation function object-template CCFdt, and 
d) Rotational profile calculated subtracting CCF2tt * G from 
CCFdt- In both panels arbitrary vertical shifts have been ap- 
plied for clarity. 

A In A m v/c, is a better approximation than the classical AA/A =s 
v/c. 

The FT of G(x) presents zeros at different positions cr„, 
which are related to v sin i by 



vsin; = 



k„(s) 



(4) 



where k n (e) are function s of the limb darkening coefficient 
(IReiners & Schmitt1l2002l) . For the first zero, k\ is given implic- 
itly by the following expression: 



- = 1 + 

S 



sin/3 cos/3 

T 2 T 



where Ji is the first-order Bessel function and fi = 2nk\ . In our 
procedure we used the following approximate formula 



0.60975 



0.0639s + 0.0205s 2 + 0.021s 3 , 



which is accurate enough (better than 0.01%) for the whole 
range of limb darkening coefficient of normal main-sequence 
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Fig. 2. Rotational profile recovered from the CCF. Upper panel: 
the functions involved in equation[2] from top to bottom: CCFdf , 
CCF27-7-*Gi, and the recovered rotational profile G. Lower 
panel: comparison of the retrieved rotational profile (filled cir- 
cles) with the input theoretical profile (solid line). 



stars: e = 0.00-1.10. We note that this cubic is not a series ex- 
pansion around s = but a fit of the function in this interval. In 
fact, it i s more accurate than th e fourth-degre e polynomial calcu- 
lated by iDravins et al. I (fl990l eq. 6; eq. 8 in lReiners & Schmitt I 
120021) in the range of stellar limb darkening coefficients. 

Once k\ (e) has been calculated, the determination of v sin / 
results from eq.|4] We calculate the power spectrum of the ex- 
tracted rotational broadening function to evaluate the position 
of the first zero en of the FT and calculate vsin/ from eq.|4] 
As mentioned before, this procedure is valid for stars rotating 
rigidly. The eventual presence of differential rotation might be 
detected by studying the shape of th e derived rotational profil e 
or using the first two zeros of the FT dReiner s & Schmitt 2002). 

Finally, it is important to mention that a linear limb dark- 
ening law was used in this work, but it is possible to consider 
other types of limb da r kenin g law like the ones compared by 
Bro wn & Verschueren I (1 19971) . We estimate that a linear law 
could produce errors of ~ 1%; therefore, a more adequate limb 
darkening law is a potential improvement for the future. 



3. Procedure and assumptions 

The procedure was programmed as an IRAF task divided in five 
stages: 



1. calculation of the CCF between an object spectrum and a 
template spectrum, 

2. extraction and cleaning of the CCF central maximum, 

3. determination of the first zero of the FT of the CCF maxi- 
mum, 

4. calculation of v sin i, 

5. measurement error estimation. 

In this section we describe four important aspects for implement- 
ing our method. 

3.1. Influence of nonrotational broadening effects 

We evaluated the impact of nonrotational broadening effects 
on the position of the first zero and on the shape of the main 
lobe of the FT using a synthetic spectrum for a J e ff = 9000 
K atmosphere model in the wavelength range 4480 - 4600 
A. This spectrum was convolved with a rotational profile of 
vsin/ = 40 kms -1 and different Voigt profiles to simulate ad- 
ditional nonrotational broadening. The FT of the various broad- 
ening profiles, derived using the CCF, are compared in Figure[3] 
We used Voigt profiles corresponding to different combinations 
of Gaussian (gfwhm) and Lorentz (lfwhm) components' full- 
width-half-maximum(FWHM). The change in the general shape 
of the first lobe with the nonrotational broadening effects is more 
significant than the shift of the first root, which in the four ex- 
amples lies within 0.4 km s _1 fr om the original rotat i onal p rofile. 
This has been already noted bv lReiners~ & Schmitt] d2002l) . 

We estimate that acceptable vsin / measurements are ob- 
tained as long as the rotational profile is at least twice wider than 
the nonrotational profile. In main-sequence stars, typically the 
intrinsic FWHM of metallic lines is 0.07-0.12 A. Therefore, for 
high-resolution spectra, the lower limit for v sin i measurements 
through the zero of the FT is about vsin; =5-8 kms -1 . B-type 
stars are more difficult to measure since their most conspicuous 
spectral features are He i lines, whose intrinsic widths are usu- 
ally on the order of 0.8 A. Such lines would not be suitable for 
rotation measurement with this technique, except for very fast 
rotators (vsin / > 100km s _1 ). In fact, these lower limits depend 
on the S/N, since the effect of the intrinsic line broadening is to 
reduce the intensity of the first subsidiary lobe, making it more 
difficult to determine the first zero position. 

All in all, for the usual values of S/N and instrumental broad- 
ening, the variation in the first zero position caused by additional 
broadening and noise is below 1%. A typical case is shown in 
Fig|4] where the FT of 20 spectra with different random noise 
corresponding to S/N = 100 are plotted. The atmospheric pa- 
rameters, vsin/, and spectral range are the same as for Fig. [3] 
and the assumed instrumental profile is 6.6 kms -1 . We mea- 
sured (vsin/) = 39.88kms-'with o- = 0.26 kms -1 (0.7%) for 
S/N = 100. Using S/N = 50 the typical error is cr = 1.3%. 
Considering that we used just a 120 A wide region, we con- 
clude that, with medium quality spectra (e.g. A/l « 1000 A and 
S /N - 50 or AA ss 100 - 200 A and S/N = 100), errors are well 
below 1%. 

In light of these results and under the assumption of rigid ro- 
tation, we consider that a global fitting of the FT of the rotational 
profile does not present much of an advantage for measuring ro- 
tational velocity, since the general shape of the function depends 
on the nonrotational broadening effects much more strongly than 
the position of the first root does, and therefore, additional free 
parameters should be included in the fitting in order to model 
the nonrotational contributions simultaneously. Moreover, even 



4 



C. G. Diaz et al.: Accurate stellar rotational velocities based on CCFs 



vsini (km s ) 
' 50 40 30 ao 



42 40 38 





III 1 


1 1 


iii i - 










>, 






















log Int 






















; 






-4 




i \i 







10 20 30 40 16 17 

a x 1000 (km"' s) 

Fig. 3. Fourier transform of the CCF of 5 synthetic spectra 
(AA = 120A) with r eff = 9000 K, vsin i = 40 kms -1 and con- 
volved with the following Voigt profiles: Thick solid line: orig- 
inal rotational profile. Thin solid line: gfwhm = 0.2A. Dotted 
line: gfwhm = 0.1A and lfwhm = 0.1A. Dashed line: gfwhm = 
0.4A. Dot-dashed line: gfwhm = 0.2A and lfwhm = 0.2A. 
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Fig. 4. Fourier transform of the CCF of 20 synthetic spectra 
(AA = 120A) with T eff = 9000 K, vsin / = 40 kms" 1 and ran- 
dom noise simulating S/N = 100. 



though the S/N is evidently lower around the first root, it is usu- 
ally high enough for this purpose, since the CCF peak has a 
much higher S/N than individual spectral lines. 

We note that the FT zero position is even less sensitive 
to nonrotation al broadening than the Bess e l-Fou rier transform 
proposed by iPiters. Groot. & van Paradijs I (119961) . In fact, ac- 
cording to Fig. 2 of Piters et al. (1996), the maximum of the 
Bessel-Fourier transform is shifted by about 1% (2%) when 
an additional Gaussian broadening of gfwhm = 0.14 x vsin/ 
(0.28 x v sin /) is present. These values approximately correspond 
to 0.2 and 0.4 A in the calculations of our Fig. [3] for which we 
found a shift in the FT zero of only 0.2% (0.4%). 



3.2. Considerations on limb darkening 

The influence of the limb darkening coeffici ent on the rotational 
broadening function G(A) is we ll known (see lReiners & Schmitt I 
l2002t ICollins & Truax 1 1 19951) . Although the errors can rise 
to n early 16 %, most autho r s assume a fixed value of s = 
0.6 (ISimon-Diaz &Herrerol l2007t ISimon-Diaz et al. I 12 006: 
Reiners & Rover 2004; Rover et al. I l2002allbt Collins & Truax 
119951: iRamella et al. 1 1 1989b . even for low-mass stars (spec- 
tral types M or L) and brown dwarfs (IWhite & Basri 1 12 003; 
i Mohantv & Basri 1 120031: iBailer-Jonesl l2004t iTinnev & Reid I 
1998). 

It is also worth noting that the limb darkening coefficient in 
the core of the lin es may differ substantia lly from that of the 
nearby continuum. ICollins & Truaxl (1 19951) treated this problem 
and indicate that the difference may be important. For a non- 
rotating B9V model they calculated e = 0.57 for the contin- 
uum at A — 4475A and s = 0.29 for the core of the Mg n line 
at /14481.13. Even though the limb darkening coefficient varies 
within the intrinsic line profile, it is valid to define the rotational 
profile G(x) as in eq. [3] using an effective limb darkening coeffi- 
cient (e.g. weighted average over the line profile), provided the 
rotational broadening is much larger than the intrinsic line pro- 
file. Our method does not use single lines but spectral regions 
that include absorption lines with s values that differ from the 
values of s for the continuum in different amounts. This can still 
be considered by defining an effective limb darkening coefficient 
for each spectral range. However, we do not consider this issue 
in the present paper. Instead, we split the spectrum to perform 
the measurements in regions of about 200-400 A and adopt the 
limb darkening coefficient for the continuum at the central wave- 
length of each region as representative of the spectral region. 

Moreover, owing to the dependence of e on temperature and 
wavelength, the parameter k„(s) in eq.|4]varies with spectral type 
and along the spectrum. To account for this effect, we include 
a calibration k\ (A, T e ff ) (see next section) specifically obtained 
and tested for A type stars. The purpose of having an empirical 
calibration for a particular range of temperatures or spectral type 
is to maximize the precision and to keep the procedure as simple 
as possible. Another effect that has not been considered is the 
gravity darkening, which is important for v sin / larger than 100 
kms" 1 . ICollins & Truax I (1 19951) also discusses this effect and 
warns that the error may be as large as 10% for v sin i larger than 
200 kms -1 . Consideration of this effect may also be a second 
step towards improving the method. 

3.3. Template spectrum and spectral region selection 

The template spectrum must be morphologically similar to the 
object spectrum but with zero rotational velocity. As the maxi- 
mum of the CCF is the result of the contribution of all the co- 
incident spectral lines between both spectra, slight differences 
resulting from small discrepancies in spectral type, metallicity, 
or the presence of spectral peculiarities have little influence on 
the CCF maximum. We tested the incidence of using templates 
that differ morphologically from the target spectrum, measur- 
ing a synthetic spectrum for an A7V star (r e ff = 8000 K and 
logg=4.0) and vsin/ = 60 km s~' with various templates. We 
measured six spectral regions > 240 A in the range 3985 - 5870 
A with 19 main sequence templates ranging from r e jf =5 000 
to 14000 K, and we found that the error produced by spectral 
type mismatch is below 1 % for templates in the range AO - GO 
(r eff =6 000 - 10 000 K) as shown in Figure Although the 
present work is focused on A-type stars, from Fig.|5]it is also ev- 
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Fig. 5. Error in (v sin i) from the mismatch of spectral type be- 
tween the object and the template spectrum. The .x-axis has the 
temperature of the template. The temperature of the object is in- 
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Fig. 6. Percentage error on (v sin i) from single regions of arti- 
ficially broadened spectra. For each velocity two values of S/N 
are plotted. S/N = 100 and 200 for (vsinz) = 80, 180, and 280 
km s- 1 . S/N = 100 and 70 for (v sin z) = 20, 120, and 220 km s" 1 . 
Triangles: Standard deviation of 10 consecutive measures (true 
error). Circles: Designated error (calibration). 



ident that earlier spectral types are more sensitive to this effect, 
so this analysis should be carried out for the spectral range of 
interest in future works. 

We also used an A5 V template to measure AO and A9 stars in 
our program and a dwarf template to measure giants and super- 
giants of the same type. In all cases, the differences in the v sin i 
are always below 1%. Such a small error shows the strength of 
the CCF, and that spectral type mismatch is not the main source 
of error in A-type stars. For this error to be comparable to the 
dispersion of values obtained from different regions in the same 
spectrum, the error in the template temperature has to be > 2 000 
K. 

The selection of the spectral regions to be correlated depends 
on the quality and the wavelength range of the observed spectra. 
Since the central wavelength of each region is used to calculate 
k\ , we consider 500 A as the maximum region width to minimize 
the influence of the limb darkening coefficient on k\ . Typically, 
in an interval of 500 A, the limb darkening coefficient varies 
within +0.015, which affects the measured rotational velocity in 
about ±0.25%. In practice, the resulting error would be signifi- 
cantly lower since it depends on the difference between coeffi- 
cient for central wavelength (used for the calculations) and the 
effective coefficient, which depends on the spectral line distribu- 
tion in the region under consideration. 

In addition, since the method is based on the rotational ve- 
locity as the main broadening factor in all spectral lines, any 
line strongly deviated from a rotational profile must be avoided. 
Otherwise, it will introduce an unwanted distortion into the CCF. 
Even though some of them might be known a priori, e.g., Hei 
in early B stars or Ca n in late A to F stars, it is convenient to 
evaluate whether other regions introduce a significant distortion 
into the CCF. The cross-correlation between two template spec- 
tra in different spectral regions can be used to detect regions that 
introduce significant distortions in the CCF. 

3.4. Single measurement error calculation 

The error assigned to a single measurement is calculated from 
the intensity and the FWHM of the CCF central maximum, and 
the noise present in the CCF. To calibrate the error formulae we 



used synthetic spectra artificially broadened with six different 
rotational profiles and random noise simulating S/N in the range 
S/N = 70 - 200. Then, ten spectra were generated for each pair 
of parameters and four regions were measured in each spectrum. 
For each region the standard deviation of the ten measures was 
adopted as the true error. 

Finally, we calibrated the error as a function of the inten- 
sity I (height of the CCF peak), the FWHM, and the noise 
of the CCF (rms). As a result, we obtained the following ex- 
pression to calculate the error of a single region measurement: 
Error(km/s) = 4.42 FWHM 520 ■ rms ■ I" 108 , valid for spectra 
with S/N = 70-200, the range used for the calibration. Figure [6] 
shows no systematic difference among real error and its calibra- 
tion. 

4. Application 

4.1. Application to A-type stars 

As part of our current research on rotational velocity, we have 
been engaged for several years in a program to measure v sin i 
for all southern A-type stars of the Bright Star Catalogue (BSC). 
This sample includes more than 800 stars . 

To test our method, we measured (vsinz) for 251 stars in 
the spectral range Al-AfQ, the results of which are presented 
in Table ?? . Amon g them 155 stars were also measured by 
iRover et al. I d2002ah . In a forthcoming paper we will publish the 
complete results of v sin z measurements for the whole sample of 
southern A-type stars of the BSC. 

The southern A-type stars of the BSC were observed spec- 
troscopically using the bench echelle spectrograph (EBASIM) 
fed with fibers coming from the Cassegrain focus of the 2.1m 
telescope at CASLEO. The spectrograph has been described by 
iPintado & Adelmanl (120031) . The resolving power of the spectra 
varies from ~ 30 000 in the blue end to about 25 000 in the red 
end. 

The observing material was reduced using IRAF packages. 
Templates of spectral type A1-A5 with resolving power 500,000 

1 The spectral types were taken from the SIMBAD astronomical 
database: http://simbad.u-strasbg.fr/simbad 
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Fig. 7. Left: Percentage error of the 251 measured stars. 
Asterisks: £i (%), individual region error indicator. Triangles: 
^2 (%X measurements dispersion indicator. Lines represent con- 
stant values of £, in kms -1 . Solid: 1 kms -1 . Dashed: 2 kms~'. 
Dot-dashed: 5 kms -1 . Dotted: 10 kms . Right: Histogram of 
the percentage error. Solid line: £ i (%). Dotted line: £2 (%) 
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Fig. 8. Position of the first zero of the FT of the rotational pro- 
files for 8 spectral regions of 8 spectra of the star HD77370. 
Dotted lines are theoretical lines corresponding to different val- 
ues of v sin i a nd have been calculated from the limb darkening 
coefficients in lDfaz-Cordoves et al. 1 (119951) . 



were selected from the Kurucz online database (lKuruczlll99l 
to measure vsin/ in objects with luminosity class V-IV, inde- 
pendently of the presence of spectroscopic peculiarities. For Am 
stars we used the template corresponding to the spectral type of 
the metallic spectrum. 

We divided our spectral range into five regions of different 
sizes according to the density of spectral lines. Near the blue 
end of the spectra, regions of 250 A were sufficient, while in the 
red end 400 A were necessary to obtain a well-defined central 
maximum in the CCF We avoided a rather small region between 
4690 and 4730 A in all spectral types, because it introduces a 
triangular base in the CCF maximum, which makes it impossible 
to unambiguously identify the first zero in Fourier space. For the 
same reason all lines from the Balmer series were also excluded. 

To simplify the calculation of the appropriate limb darkening 
coefficient for each star and wavelength region, we implemented 
a calibration k\ (A, T e g), where A is the central wavelength of the 
spectral region. We u sed t he limb darkening coefficien ts tabu- 
lated bv lClaretl d2000l) and lDfaz-Cordoves et al. 1 dl995h for var- 
ious photometric bands. Then, we fitted s(T e s) for our tempera- 
ture range of interest, i.e. T e ff = 7500 K-l 1 500 K (spectral types 
B9-A9), for log g — 4.0 and for the filters B, V, v, b, and y, whose 
central wavelength are within the range of our spectra (4000 A - 
6000 A). The residuals of the fit for each filter lie within ±0.012 
(rms=0.006), which represents an error in v sin i ~ 0.15 %. The 
final wavelength calibration was made as a function of the cen- 
tral wavelength of each filter, giving e(T e s, A) as a result. By 
means of this calibration it is possible to calculate the k\ value 
appropriate for the star's 7^ and the central wavelength of spec- 
tral region. 

The number of regions used to determine v sin i depends on 
the quality of the spectra and the wavelength range, and is totally 
independent of v sin i. The observational material analyzed here 
consist of 32 spectra in the range 4000-6000 A in which we 
measured five regions, 116 spectra in the range 4000-5500 A in 
which four regions were measured, 99 spectra in the range 3850- 
5000 A or in the range 4000-5250 A in which three regions were 



measured, and finally four objects in which only two regions 
were measured owing to excessive noise in part of the spectra. 

Once v sin i and its error were obtained in all the regions of 
the same spectrum, the weighted mean value of v sin i was calcu- 
lated, along with two different estimates o f its uncertainty, $\ and 
£2, fo llowing the formulae proposed by iGonzalez & Lapasset I 
(2000, eq. 1) for radial velocities. 

The error £1 is the standard deviation of the weighted mean 
of a sample of n uncorrelated observations with standard de- 
r -9i-°- 5 

viations e,: %\ = £<? ( . On the other hand, £2 is a gen- 
eralization of the standard error of the mean crj -\/n for a set 
of measurements weighted according to their individual errors. 
Therefore, £2 is calculated from the dispersion of vsin/ from 
different spectral regions, while %\ is computed from the single 
measurement errors estimated through the calibration described 
in Sect. 13.41 In general, %\ would be more appropriate when the 
number of measurements is small (2-3 spectral regions), and 
consequently the dispersion of measurements is a less precise 
estimate of the true error. In Figure [7] we have plotted the val- 
ues of £1 and £2 computed for each object. Excluding slow rota- 
tors (v sin i below 30 km s _1 ) whose line widths are close to the 
limit imposed by the spectral resolution (c/R = 10-12 kms -1 ), 
the error £1 has an average value of 1.1%, always smaller than 
2.5%. Regarding the dispersion of the values obtained from dif- 
ferent regions of the same spectrum, only three objects with 
(vsin/) > 30 kms -1 were measured with £2 over 5%, and the 
average error £2 of these objects is 1.5%. 

4.2. Application to late type stars 

One field of application for the proposed method is the measure- 
ment of rotation in stars of spectral type F or later, where it is 
not possible to measure individual lines even for low rotators. 
As a test, we measured the rotational velocity of HD 77370, an 
F3 V star with projected rotational velocity of about 60kms~ 1 , 
high enough for all spectral lines to appear blended. Eight high- 
resolution HARPS spectra of this star were downloaded from 
the ESO archive database. These spectra cover a spectral range 
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3800-6900 A with a S/N * 200. Rotational velocities were 
measured for eight spectral regions of about 200^1-00 A, using as 
template a synthetic spectrum with T e ff=6750K and solar chem- 
ical abundances. Figure[8] shows the measured vsin/ values for 
the eight spectra. 

The dependence of the position of the first zero of the FT 
of the rotational profile with wavelength - due to the variation 
of the limb darkening coefficient with wavelength - is notice- 
able. The rotational velocity of each spectrum was calculated as 
the average of the eight spectral regions. The mean value of the 
64 measurements is 59. 16km s~'. The standard deviation of the 
measurements of the 8 spectra in one single spectral region is on 
average 0.41 kms -1 . If we take this value as the typical uncer- 
tainty of one of the 64 individual measurement, the uncertainty 
of the mean v sin i for each region should be V8 times lower 
(0.14kms~'). However, if we calculate the average for each of 
the eight spectral regions, then the dispersion of these mean ro- 
tational velocities is 0.37 km s~', indicating that the main source 
of uncertainty in these measurements is not random, but small 
systematic differences between the spectral regions. All in all, 
the uncertainty of the v sin i measured from a high-quality spec- 
trum of an F-type star would be about 0.6%, demonstrating that 
high-precision rotational velocities can be measured even when 
lines are blended. The key point in this respect is the use of the 
template-template correlation function to model the subsidiary 
lobes of the object-template CCF during the calculation of the 
rotational profile. As an illustration, Figure [9] shows the recon- 
struction of the rotational profile for a small spectral region (70 
A) of the spectrum of HD 77370, where all spectral lines are 
blended. 




5. Discussion 

To evaluate the performance of our metho d in determining v sin ;', 
we compared our results with those from lRover et al. I (|2()02a). 
They used FT of line profiles to provide accurate v sin i of a large 
sample of A-type stars in the southern hemisphere observed with 
similar resolving power (~ 28000). FigurefTOlshows Royer et al. 
values of v sin i and the ones obtained in this work for 155 objects 
in common with their sample. Even though a good agreement 
is found in stars of moderate, projected rotational velocities, a 
significant deviation is noticed for vsinz' > 150 km s~' with an 
average difference of 5 ± 1%, our values larger than theirs. 

There are several dif ferences between the methodology used 
bv iRover et al. I (1200 2a) and the present work that can account 
for the effect present in Figure [TOj First, they measured single 
lines in the range from 4200 A to 4500 A, i.e. a region of 300 A 
wide which is the mean width of each region used in this work. 
Second, the blending effect becomes evident through the depen- 
dence of the number of measured lines with v sin i. They used 
in average less than 3 lines on objects with vsin/ > 60 kms -1 , 
and over 90 km s~'the average number of measured lines is < 2. 
Even though the dependence on the position of the first zero of 
the Fourier transform with the limb darkening coefficient was not 
taken into account by Royer et al., the adopted value e = 0.6 is 
typical of the limb darkening coefficient in the range of temper- 
atures of the objects in the sample. However, an error of ±0.01 
in k\ represents an error of 1.5% in vsin/. Thus, using a wrong 
k value could account for 1.5% of the deviation detected in the 
comparison. 

Finally, the mai n source of e r ror is t he uncertainty on the 
continuum position. IRover et al. I d2002al) analyzed the magni- 
tude of this effects by broadening synthetic spectra of T e ff from 



Fig. 9. Example of the rotational profile calculated from a high- 
S/N spectrum of only 70 A for an F-type star. Upper panel: 
Object spectrum and template spectrum. Lower panel: the same 
as Fig. 1. 
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Fig. 10. Comparison between our results and IRover et al. I 
(2002a) for objects in common. Squares: normal stars. Circles: 
peculiar stars. Solid line: Normal stars fit /((vsini)) = 
1.43 (vsin i) 092 



7500 to 10000 K with vsin; = 10, 50, and 100 kms 1 . They 
found that in the coldest models with the largest v sin i the error 
can reach 3%. This uncertainty affects spectral lines giving as 
a result a lower value of v sin i than the real value. To evaluate 
the intensity of this type of systematic error with our method, 
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Fig. 11. Difference between real vsini and average over ten 
consecutive measures with distinct random noise. Dotted line: 
/«v sini» = 0.0032 <v sini) + 0.0456. 



we measured vsini in a synthetic A5V spectrum broadened by 
vsim = 20, 80, 120, 180, 220, and 280 kms -1 . Random noise 
was added to simulate S/N = 100, generating ten spectra with 
different noise pattern for each value of v sin i. Then, four regions 
were used to measure v sin i, and the average over the ten mea- 
sures for each region was calculated. As expected, we found a 
tendency to measure a lower value of v sin i than the actual value. 
This effect is shown in Figure[TT] Nevertheless, with the method 
proposed in this work the difference between measures and real 
values stays under 1% (( "T^O iso-aootms-- = 0.997+0.004) 
even for vsin/ = 280 kms -1 . In conclusion, a considerable 
fraction of the diffe rence between our results and those from 
Rove r et aT71 (|2002a) can be attributed to the effect of line blend- 
ing on the continuum position. The influence of the line blending 
in high rotational velocity stars is significantly less with the pro- 
posed method than with Fourier transform of single-line profiles, 
by virtue of subtracting the secondary maximums from the CCF 
before calculating the Fourier transform. Similar results would 
be obtained by removing the blended lines with a similar proce- 
dure in the spectrum: modeling the smaller lines by convolving 
a rotational profile with a synthetic spectrum in which the line 
of interest has been subtracted. 

The main originality of the proposed method is the way in 
which the broadening function is built from the CCF, which, on 
the one hand, allows measurement of stars with blended spectral 
lines and, on the other, improves the S/N. Once the broaden- 
ing function has been obtained, we determine the stellar rotation 
through the first zero of the FT. However, other alternative analy- 
ses are possible, such as direct fitti ng of the broadening function 
or Fo urier-Bessel transformation ( Piters. Groot. & van Paradrjsl 
Il996h . 

As an illustration, we used our method in combination with 
the Fourier-Bessel transformation to measure one of the objects 
in Piters et al. (1996) for which a HARPS spectrum was publicly 
available. These authors applied the Fourier-Bessel transforma- 
tion method to individual spectral lines for measuring rotation 
in F-type stars. Owing to line blending, stars with vsin; > 100 
kms 1 could not be measured or had errors of about 10-20 
%. We selected one of their fast rotating objects: HD 1231 1, for 
which they measured 110+22 kms -1 . We defined 5 spectral re- 
gion between 3980 A and 6530 A obtaining an average of 157.4 
km s 1 with cr = 4.6 km sousing our method with the FT first 



zero, and 155.9 kms~' with cr = 1.8 km s~' from our broadening 
function but using the Fourier-Bessel transform. Therefore, the 
application of Fourier-Bessel transform method to the broaden- 
ing function derived from the CCF would be an excellent alterna- 
tive for late-type fast-rotating stars. In the case of slowly rotating 
stars, however, the Fourier-Bessel transform is more sensitive to 
broadening effects that are different from rotation, so the zero of 
the FT would give better results (see sect. 3.1). 

6. Summary 

We have developed a method for measuring v sin i based on the 
FT of the central maximum of the CCF. This combination pro- 
vides a simple and precise solution to the line blending problem 
with medium-resolution spectra and/or spectral types later than 
mid-A. As a result, the number of useful spectral regions per ob- 
ject with our method is independent of spectral type and v sin i. 

The high precision of the proposed method is supported by 
two key features of the procedure. The first one is to use an em- 
pirical calibration to take limb darkening effects into account in 
the position of the first zero of the FT, which is directly related to 
v sin i. The second key feature is the subtraction of the subsidiary 
lobes of the CCF during the calculation of the rotational broad- 
ening function. This assures a good definition of the rotational 
profile even in presence of line blending. 

The v sin i value is derived from the first zero of the FT of 
broadening function. However, the S/N of the resulting rotational 
profile is high enough to also do other detailed studies of the 
shape of the reconstructed rotational profile. 

For slowly rotating stars, a lower limit for measurable v sin i 
is imposed by the instrumental profile and the intrinsic line 
width, which is included twice in our broadening function owing 
to the use of a template with lines of finite width. For metallic 
lines this limit is 5-8 km s . In early-type stars the measurement 
is more difficult since metallic lines become weaker and He lines 
are usually not useful. Morphological differences between tem- 
plate and object are also less significant in late-type stars. For 
late-A type stars or later, differences in several subtypes have no 
impact on the rotation measurements. 

We consider that, in accurate v sin i determinations, one sig- 
nificant contribution to the errors is the adopted value for the 
limb darkening. Even though in this work we use a wavelength- 
dependent coefficient, the continuum limb darkening might dif- 
fer significantly from the line limb darkening. Using many spec- 
tral lines in our calculations, this difference might be averaged 
out, but systematic differences between lines and continuum 
might result in small but systematic v sin i errors. This is an issue 
that deserves further research. 

We applied the proposed method to a sample of 25 1 A-type 
stars. Measurement errors are always under 2.5%, being 1.1% 
the average error for stars with (vsin/) over 30 kms 1 . As re- 
gards the dispersion of values measured on different regions 
across a spectral range of ~ 11 50 A to 2000 A wide, we found 
standard errors of the mean below 5% with an average of 1.5%. 
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Table 1. Results for 251 A1-A5 stars from the Bright Star Catalogue in the southern hemisphere. 
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Table 1. Continued. 



HD 


(vsini) 




e 2 


Spectral 




km .s ,_1 


km s l 


km s' 1 


Type 


66210 


30.9 


1.5 


0.5 


A2V 


68862 


91.9 


0.5 


1.0 


A3V 


69665 


98.2 


0.9 


2.6 


A1V 


70340 


10.6 


0.3 


0.2 


A2VDn-EuSrCrSi 


70612 


104.1 


0.6 


0.7 


A3V 


71267 


19.8 


0.3 


0.4 


A3m 


71688 


132.1 


1.6 


0.7 


A1/A2V 


71815 


39.1 


0.3 


0.6 


A1/A2V 


72660 


8.1 


0.4 


0.2 


A1V 


72968 


15.9 


0.5 


0.4 


AlpSrCrEu 


73997 


193.7 


2.8 


2.2 


AlVn 


74190 


59.4 


0.4 


0.7 


A5m 


74341 


97.6 


0.9 


1.5 


A3V 


74879 


63.5 


0.6 


0.8 


A3IV/V 


75171 


96.5 


0.6 


0.6 


A4V 


75630 


163.7 


1.6 


2.1 


A2/A3IV 


75737 


16.5 


0.4 


0.7 


A4m 


75926 


194.5 


2.2 


1.4 


AlVn 


76483 


68.9 


0.4 


0.5 


A3IV 


78045 


30.6 


0.4 


0.5 


kA3hA5mA5v 


78676 


56.1 


0.4 


0.5 


A4IV 


78922 


108.3 


0.5 


0.3 


A4IV-V 


79193 


10.4 


0.5 


0.2 


A3III:m+A0V: 


80007 


145.7 


1.2 


2.2 


A2IV 


80447 


15.1 


0.2 


0.3 


A2Vs 


80951 


23.7 


0.9 


0.1 


A1V 


81157 


36.6 


0.6 


0.7 


A3IVs... 


81309 


10.6 


0.2 


0.7 


A2m 


82068 


258.9 


3.5 


2.8 


A3Vn 


82446 


54.0 


0.6 


1.6 


A3V 


83520 


149.6 


1.3 


1.8 


A2,5V 


83523 


134.4 


0.9 


1.8 


A2V 


85558 


134.6 


1.4 


3.4 


A2V 


86266 


186.6 


1.6 


0.9 


A4V 


86301 


218.9 


1.7 


2.0 


A4V 


88024 


92.3 


1.3 


3.1 


A2V 


88372 


259.2 


4.2 


8.3 


A2Vn 


88522 


23.1 


0.5 


0.4 


A1V 


88842 


78.9 


0.6 


0.3 


A3IV-V 


88955 


103.6 


0.9 


1.1 


A2Va 


89263 


67.2 


0.8 


0.5 


A5V 


89816 


214.3 


2.0 


6.6 


A4IV/V 


90630 


90.9 


0.7 


0.4 


A2,5V 


90874 


65.7 


0.5 


0.4 


A2V 


91790 


111.7 


0.8 


1.3 


A5IV/V 


93397 


94.0 


0.4 


0.9 


A3V 


93742 


49.8 


0.5 


0.9 


A2IV 


93903 


19.4 


0.4 


0.6 


A3m 


93905 


118.3 


1.2 


2.2 


A1V 


94985 


181.0 


1.4 


4.1 


A4V 


95370 


111.6 


1.1 


2.5 


A3IV 


96124 


156.3 


1.9 


0.8 


A1V 


96146 


7.7 


0.3 


0.03 


A2V 


96441 


130.3 


1.0 


1.3 


A1V 


96723 


25.4 


0.5 


0.1 


A1V 


103101 


42.1 


0.5 


0.4 


A2V 


103266 


165.1 


1.6 


1.0 


A2V 


104039 


12.6 


0.3 


0.2 


A1IV/V 


105776 


153.1 


1.1 


5.0 


A5V 


105850 


126.8 


1.2 


1.4 


A1V 


106819 


79.0 


0.6 


1.1 


A2V 


107070 


202.3 


1.9 


3.0 


A5Vn 


108107 


242.1 


3.2 


5.2 


A1V 


108925 


188.9 


1.7 


2.7 


A3V 


109074 


89.9 


0.6 


1.2 


A3V 


109704 


155.4 


1.3 


1.8 


A3V 
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Table 1. Continued. 



HD 


(vsini) 


E i 


e 2 


Spectral 




km s~ l 


km s~ i 


km s^ 1 


TvDe 


109787 


296.8 


3.8 


7.3 


A2V 


111588 


60.6 


0.5 


0.5 


A5V 


114330 


7.4 


0.5 


0.2 


AlIVs+... 


114576 


162.0 


2.0 


4.2 


A5V 


115892 


90.3 


1.0 


1.6 


kA15hA3mA3va 


116061 


194.2 


1.5 


2.4 


A2/A3V 


116197 


239.7 


3.3 


4.4 


A4V 


117150 


220.4 


2.4 


5.9 


A1V 


117558 


169.8 


2.2 


2.1 


A1V 


118098 


232.6 


2.2 


6.3 


A3V 


119938 


69.7 


0.4 


0.3 


A3m... 


122958 


172.2 


1.9 


7.1 


A1/A2V 


123998 


17.2 


0.3 


0.1 


A2m 


124576 


126.5 


1.4 


2.1 


A1V 


125283 


257.3 


3.6 


4.0 


A2Vn 


126367 


65.8 


0.6 


0.7 


A1/A2V 


126722 


101.2 


0.7 


1.1 


A2IV 


127716 


96.4 


0.8 


0.3 


A2IV 


130841 


59.6 


0.9 


0.8 


kA2hA5mA4Iv-v 


132219 


42.9 


0.4 


1.8 


A4V+„, 


133112 


96.4 


0.6 


1.0 


A5m 


134482 


160.4 


1.0 


1.5 


A3IV 


134967 


259.4 


3.2 


13.2 


A2V 


135235 


75.9 


0.5 


0.3 


A3m 


135379 


68.5 


0.9 


0.8 


A3Va 


137015 


119.6 


1.4 


3.4 


A2V 


137333 


44.1 


0.6 


0.8 


A2V 


138413 


21.4 


0.2 


0.1 


A2IV 


138965 


102.7 


1.0 


2.0 


A1V 


141413 


116.6 


0.6 


0.8 


A5IV 


142139 


89.2 


0.7 


0.7 


A3V 


142445 


137.7 


1.1 


1.4 


A3V 


142629 


43.0 


0.6 


1.1 


A3V 


143101 


86.1 


0.6 


0.6 


A5V 


145570 


42.3 


0.6 


0.5 


A3IV 


145607 


217.7 


2.7 


4.2 


A4V 


145689 


106.4 


0.8 


0.5 


A4V 


146667 


225.2 


2.5 


1.6 


A3Vn 


148367 


19.5 


0.8 


0.8 


A3m 


150573 


235.8 


2.6 


1.6 


A4V 


150894 


136.2 


0.8 


0.9 


A3IV 


151676 


149.4 


1.2 


1.1 


A3V 


152127 


31.6 


0.5 


0.7 


A2Vs 


153053 


102.8 


0.7 


1.0 


A5IV-V 


154310 


263.0 


3.5 


7.7 


A2IV 


154418 


72.6 


0.6 


0.6 


Aim,,, 


154494 


114.8 


0.6 


0.8 


A4IV 


154895 


121.6 


1.1 


4.0 


A1V+F3V 


155259 


214.0 


3.2 


4.0 


A1V 


159492 


54.1 


0.4 


0.7 


A5IV-V 


160613 


112.6 


1.2 


1.6 


A2Va 


164577 


194.6 


1.8 


3.5 


A2Vn 


165189 


129.8 


0.8 


0.7 


A5V 


165190 


129.2 


0.7 


1.4 


A5V 


166393 


191.7 


3.0 


3.0 


A2V 


166960 


27.3 


0.3 


0.3 


A2m 


167564 


153.0 


1.3 


1.8 


A4V 


169853 


22.6 


0.3 


0.7 


A2m 


170384 


131.3 


0.6 


0.5 


A3V 


170479 


67.6 


0.7 


2.2 


A5V 


170642 


177.3 


1.3 


5.3 


A3Vn 


171856 


106.1 


0.8 


0.2 


A5IV 


172555 


116.4 


0.6 


0.5 


A7V 


172777 


148.6 


1.6 


4.6 


A2Vn 


173715 


122.3 


0.9 


4.6 


A3V 


175638 


151.1 


0.9 


2.7 


A5V 
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Table 1. Continued. 



HD 


{vsini) 




e 2 


Spectral 




km 


km 


km s ^ 


TvDe 


175639 


214.3 


1.8 


2.9 


A5Vn 


176687 


68.9 


0.4 


0.3 


A2,5Va 


176984 


30.0 


0.6 


0.1 


A1V 


178253 


203.2 


1.7 


2.8 


A2Va 


180482 


80.4 


0.7 


1.4 


A3IV 


181383 


180.7 


2.3 


1.4 


A2V 


183545 


226.5 


2.3 


3.0 


A2V 


184552 


13.1 


0.3 


0.7 


Aim 


184586 


197.7 


2.3 


1.7 


A1V 


187421 


170.0 


2.3 


2.5 


A2V 


187653 


248.2 


3.0 


6.0 


A4V 


188899 


69.7 


0.6 


1.1 


A2V 


189118 


45.2 


0.5 


0.6 


A4/A5IV 


189388 


248.4 


3.1 


3.3 


A2,5V 


189741 


137.9 


1.2 


2.3 


A1IV 


197725 


153.4 


1.5 


3.1 


A1V 


198001 


115.4 


1.2 


2.4 


A1,5V 


198160 


183.0 


2.4 


3.3 


A2,5IV-V 


198161 


176.7 


1.9 


5.5 


A2,5IV-V 


198529 


275.4 


3.1 


2.7 


A2Vn 


198743 


51.7 


1.0 


1.2 


A3m 


199475 


212.0 


3.4 


1.1 


A2V 


201906 


264.4 


6.0 


14.1 


A1V 


202730 


135.6 


1.8 


2.6 


A5Vn: 


203562 


58.0 


0.7 


1.3 


A3V 


204394 


132.7 


1.1 


2.2 


A1V 


204854 


107.1 


1.0 


1.4 


A2IVv 


205765 


172.2 


2.2 


2.6 


A2V 


205811 


102.6 


1.3 


1.3 


A2V 


207155 


136.6 


1.5 


5.0 


A2V 


208321 


215.6 


1.9 


1.1 


A2Vn 


208565 


308.8 


6.4 


8.9 


A2Vnn 


210049 


307.7 


4.7 


2.4 


Al,5IVn 


210739 


189.7 


1.6 


3.7 


A3V 


212728 


250.1 


2.9 


0.8 


A4V 


213884 


165.0 


1.2 


3.1 


A5V 


214085 


163.4 


2.0 


2.0 


A3Vn 


214150 


80.9 


0.7 


1.7 


A1V 


214484 


8.5 


0.1 


0.5 


A2Vp 


215729 


170.5 


2.1 


0.7 


A2V 


216627 


83.2 


0.6 


0.3 


A3V 


216900 


67.0 


0.4 


0.8 


A3Vs 


216956 


91.6 


0.5 


0.4 


A4V 


217498 


91.9 


1.3 


2.1 


A2V 


218108 


194.5 


1.7 


1.9 


A3,5V 


219402 


154.2 


1.8 


2.0 


A3V 


221675 


57.2 


0.5 


0.1 


A2m 


221760 


22.4 


0.5 


0.5 


A2VpSrCrEu 


222602 


234.7 


4.3 


9.8 


A3Vn 


223438 


89.0 


0.5 


0.9 


A5m 


223466 


69.5 


0.9 


0.6 


A3V 


223991 


19.1 


0.9 


0.5 


A2V+F2V 


224361 


95.6 


0.8 


1.5 


A1IV 
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